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The geodesic deviation equation, describing the relative accelerations of nearby particles, and the 
Raychaudhury equation, giving the evolution of the kinematical quantities associated with deforma- 
tions (expansion, shear and rotation) are considered in the framework of modified theories of gravity 
with an arbitrary curvature-matter coupling, by taking into account the effects of the extra force. 
As a physical application of the geodesic deviation equation the modifications of the tidal forces due 
to the supplementary curvature-matter coupling are obtained in the weak field approximation. The 
tidal motion of test particles is directly influenced not only by the gradient of the extra force, which 
is basically determined by the gradient of the Ricci scalar, but also by an explicit coupling between 
the velocity and the Riemann curvature tensor. As a specific example, the expression of the Roche 
limit (the orbital distance at which a satellite will begin to be tidally torn apart by the body it is 
orbiting) is also obtained for this class of models. 

PACS numbers: 04.50.Kd, 04.20.Jb, 04.20.Cv, 95.35. +d 



I. INTRODUCTION 

The strong observational evidence that recently our 
universe has entered an accelerated expansion phase [![, 
and the astonishing result that around 95-96% of the 
content of the Universe is in the form of dark energy 
and dark matter, respectively, with only about 4-5% be- 
ing represented by baryonic matter [2J, has shown the 
limitations of standard general relativity. Hence, despite 
the remarkable success of general relativity at the solar 
system scale, on a galactic and cosmological scale grav- 
itational theories face two fundamental problems: the 
dark matter problem, and the dark energy problem, re- 
spectively. Although in recent years many different ap- 
proaches have been proposed to explain the observational 
results of cosmology, a satisfactory model has yet to be 
obtained. 

However, a promising way to explain the observational 
data is to assume that at large scales the Einstein gravity 
model of general relativity breaks down, and a more gen- 
eral action describes the gravitational field. Theoretical 
models in which the standard Einstcin-Hilbcrt action is 
replaced by an arbitrary function of the Ricci scalar i?, 
first proposed in Q , have recently been extensively inves- 
tigated. Cosmic acceleration can be explained by f(R) 
gravity and viable cosmological models can be found 
5j. For a review of f(R) generalized gravity models see 
0]. The possibility that the galactic dynamic of mas- 
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sive test particles can be understood, without the need 
for dark matter, was also extensively considered in the 
framework of f(R) gravity 0- In the context of the Solar 
System regime, local tests in the weak field approxima- 
tion are also to be retained by assuming the chameleon 
mechanism Jj|, where the mass of the chameleon scalar 
field depends on the local background matter density, i.e., 
in regions of high matter density the chameleon scalar 
field is massive. However, it has recently been shown 
that in the hybrid metric-Palatini gravitational theory, 
which consists of the superposition of the metric Einstein- 
Hilbert Lagrangian with an f(R) term constructed a la 
Palatini [9|, that even if the scalar field is very light, 
the theory passes the Solar System observational con- 
straints. Therefore the model predicts the existence of 
a long-range scalar field, modifying the cosmological and 
galactic dynamics. 

However, most of the generalizations of standard gen- 
eral relativity concentrated only on the geometric part 
of the action, and assumed that the matter part is un- 
changed, that is, in the total Lagrangian the matter 
term was considered as a simple additive term. This 
point of view severely limits the possibilities of the 
matter-geometry interaction, restricting the degrees of 
freedom of the gravitational theories. From a physical 
point of view, generalized gravitational models, involving 
curvature-matter interactions, cannot be ruled out a pri- 
ori [T^]. In this context, a generalization of f(R) gravity 
was introduced in fill ] , and extended in (l2j , by including 
in the theory an explicit coupling of an arbitrary function 
of the curvature scalar, R, with the matter Lagrangian 
density. As a result of the coupling, the motion of the 
massive particles is non-geodesic, and an extra force, or- 
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thogonal to the four-velocity, arises. This class of models 
of modified gravity can be denoted as generalized grav- 
ity models with a linear curvature-matter coupling, and 
they have been extensively studied recently [13[. For a 
review of modified gravity models with curvature-matter 
coupling see flil ]. 

Similar couplings between gravitation and matter have 
also been considered, as possible explanations for the ac- 
celerated expansion of the universe and of the dark en- 
ergy, in [l5j . An interesting extension of standard gen- 
eral relativity has also been been proposed, namely, the 
f(R,T) modified theories of gravity, where the gravita- 
tional Lagrangian is given by an arbitrary function of the 
Ricci scalar R and of the trace of the stress-energy tensor 
T . It is interesting to note that the dependence from 
T may be induced by exotic imperfect fluids or quantum 
effects (conformal anomaly). This extended theory of 
gravity may be considered as a relativistically covariant 
model of interacting dark energy. It was further argued 
that the new matter and time dependent terms in the 
gravitational field equations play the role of an effective 
cosmological constant. 

The above models were further generalized in a recent 
proposal where the gravitational action is given by an 
arbitrary function of the Ricci scalar R and of the La- 
grangian density of the matter L m pjj , and with a gen- 
eralized scalar field and kinetic term dependences (l8j |. 
Thus, f(R, L m ) represents the natural generalization of 
the models with linear matter coupling, as well as the 
most general extension of the standard Hilbcrt action 
for the gravitational field, S = J [R + L m ] ^/^gd 4 x. In 
this class of models the energy-momentum tensor of the 
matter is generally not conserved, and the motion is non- 
geodesic. In the particular case in which the Lagrange 
function of the matter is a function of the energy density 
of the matter only, the equations of motion of test parti- 
cles can be obtained, by using a variational principle. 

It is the purpose of the present paper to investigate 
some other interesting properties of the motion of test 
particles in gravity models with an arbitrary curvature- 
matter coupling. The weak field limit of the model is 
carefully analyzed, and it is shown that in first order in 
both Ricci scalar and matter energy density one obtains 
the Newtonian Poisson equation in the presence of an ef- 
fective cosmological constant. The equation of geodesic 
deviation and the Raychaudhury equation are also formu- 
lated by explicitly including in the equations the effects 
of the curvature-matter coupling and of the extra force. 
Some of the physical implications of the geodesic devia- 
tion equation, namely, the problem of the tidal forces in 
this class of models is also considered, and the general- 
ization of the Roche limit is obtained. 

The present paper is organized as follows. The field 
equations and the equations of motion of the model are 
derived in Section [TTJ The geodesic deviation equation 
(the Jacobi equation) of the test particles, as well as the 
Raychaudhury equation of the model, are obtained in 
Section HEH Some physical applications of the geodesic 



deviation equation are discussed in Section HVl We dis- 
cuss and conclude our results in Section [V] In the present 
paper we consider a system of units with 8irG = c = 1, 
and we follow the Landau-Lifshitz conventions for 
the metric signature and the definition of the curvature 
tensor. 



II. FIELD AND MOTION EQUATIONS WITH 
AN ARBITRARY CURVATURE-MATTER 
COUPLING 

The most general action for a / (i?, L m ) type modified 
theory of gravity involving an arbitrary coupling between 
matter and curvature is given by [l2l | 



S = j f(R, L m ) y/^g d 4 x , 



(1) 



where f(R,L m ) is an arbitrary function of the Ricci 
scalar R, and of the Lagrangian density corresponding 
to matter, L m . The only requirement for the function 
/ (R,L m ) is to be an analytical function of R and L m , 
respectively, that is, it must possess a Taylor series ex- 
pansion about any point. The matter energy-momentum 
tensor T^„ is defined as 



(2) 



By assuming that the Lagrangian density L m of the 
matter depends only on the metric tensor components, 
and not on its derivatives, we obtain = L m g^ v — 
2dL m /dg^. 

Varying the action with respect to the metric tensor 
g^v we obtain the field equations of the model as 



fli (R, Lm) R[iV + P^vfR (R, L m ) — — [/ (R, L m ) 

—fh m (R, L m ) L m ] g^v = ~ /i» (^) L m ) Tfj, v , (3) 

where we have denoted fa (R, L m ) = df (R, L m ) /dR 
and f Lm (R, L m ) = df (R, L m ) /dL m , respectively, and 
we have introduced the operator P^v, defined as 



(4) 



with □ = V^V 1 . The operator P MJ , has the property 



PS = 3D. 



By contracting the field equations Eq. (j3]), we obtain 
the scalar equation 



ZUf R + f R R - 2/ = ( -T - 2L r 



h 



(5) 



where T = is the trace of the matter energy- 
momentum tensor. By eliminating the term D/r (R, L m ) 
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between Eq. <|3j> and Eq. ((5j) , we can reformulate the field 
equations as 



R[iu — A (R, L m ) g^ 



1 



f R (R, L m ) 



V m V v /h(J2, L m ) 



+^(R,L m )[T liV --Tg^),{&) 



where we have denoted 



A(R, L m ) 



and 



2f R (R, L m ) R — f(R, L m ) + f Lm (R, L m ) L r , 



$(R,L m ) = 



6f R (R,L m ) 



fL m ffi; Lm) 

fn (R, L m ) 



(7) 
(8) 



respectively. 

By taking the covariant divergence of Eq. © , we ob- 
tain for the divergence of the energy-momentum tensor 
T MI/ the following relationship 

VT^ = V M In [fh m (R, L m )] (L m g^ u — T^ u ) 

= 2V^n[/ im (i?,i m )]^. (9) 

Now, assuming that the matter Lagrangian is a func- 
tion of the rest mass density p of the matter only, from 
Eq. (j9|) we obtain explicitly the equation of motion of the 
test particles in the / (R, L m ) gravity model as 



D 2 x fl 



d 2 x» 
~ds 2 ~ 



vjj\ 



F, (io) 



where the wordline parameter s is taken as the proper 
time, = dx^ /ds is the four-velocity of the particle, 
T v a p are the Christoffel symbols associated to the metric, 
and the extra-force / M is defined as 



-V„ In 



fh m (R, L Tl 



dL m (p) 



dp 



(u^v - sH 



(ii) 

The extra- force / M , generated by the curvature-matter 
coupling, is perpendicular to the four-velocity, f^u^ = 0. 
Due to the presence of the extra- force / M , the motion of 
the test particles in modified theories of gravity with an 
arbitrary coupling between matter and curvature is non- 
geodesic. From the relation U^V V U^ = it follows that 
the force / p is always perpendicular to the velocity, so 
that U^p = 0. 

The generalized equations of motion Eq. (1101) can be 
derived from the Kahil-Bazanski Lagrangian [20J, 



L p = 9^^£- + f,if = g^Wrf + frf, (12) 

where rf is an arbitrary four-vector, which can be taken, 
for example, as the deviation vector (see Section ITV|) . and 
we have denoted, for simplicity, 



ds 



drf 
ds 



(13) 



Then we obtain immediately dL p /diY = U a , and 
dL p /d?f = r u al3 U u UP + f a , respectively. Finally, the 
Lagrange equations 



d_ ( dLp\ _ dip _ 
ds \dfj" J dtf ~ ' 

provide the equations of motion Eq. (|10l) . 



(14) 



III. WEAK FIELD LIMIT OF THE FIELD 
EQUATIONS IN f(R,L m ) GRAVITY 

Generally the matter Lagrangian L m is a function 
of the matter energy density p, the pressure p as well 
as the other thermodynamic quantities, such as the spe- 
cific entropy s or the baryon number n, so that L m = 
L m (p,p 7 s 7 n). In the simple (but physically the most 
relevant) case in which the matter obeys a barotropic 
equation of state, so that the pressure is a function of 
the energy density of the matter only, p = p(p), the mat- 
ter Lagrangian becomes a function of the energy density 
only, and hence L m = L m (p) . Then, the matter La- 
grangian is given by 



12] 



Lrn (p) = P 1 + 



P dp 
P 



P(P)- 



(15) 



while the energy-momentum tensor can be written as 



= [p + P (p) + P n ( P )} u»u v - P ( P ) <r 



respectively, where 

n(p) = 



;dp 



p dp p(p) 
o P 



P 



(16) 



(17) 



The expression Tl(p) +p(p)/p represents the specific 
enthalpy of the fluid. From a physical point of view II (p) 
can be interpreted as the elastic (deformation) potential 
energy of the body, and therefore Eq. (|16j) corresponds 
to the energy-momentum tensor of a compressible elastic 
isotropic system. 

Next, we consider the weak field limit of the gravita- 
tional field equations. First, wc assume that p 3> p, and 
therefore we systematically neglect the pressure term. 
Then, from Eq. (fT"5"|) it follows that L m = p, so that 
the energy- momentum tensor is given by = pUfJJu. 

Secondly, we consider non-relativistic macroscopic mo- 
tion. Consequently we can neglect the spatial compo- 
nents in the four- velocity, and retain only the time com- 
ponent, so that U 11 = Up w (1, 0, 0, 0). Since in the weak 
field limit one can omit all time derivatives as well as the 
terms containing the products of the Christoffel symbols, 
we obtain for i? o t ne expression [l9j 



'o,iO = T% = d<t>/dx\ i = 1,2,3 and 

R\ = d 2 (j>/ (dx 1 ) , i = 1,2,3 (no summation upon the 
index i). Therefore we obtain R sa 2A(f>. 



Similarly, we have T, 
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Now, multiplying Eq. (|6|) with U^U 11 gives the scalar 
equation 



fa (Rj L m ) 



+A (R, L m ) + $ (R, L m ) (u»U u T^ - M . (19) 



In the weak field limit and for a static geometry, 
Eq. (|19j) immediately provides the following relationship 

A<j> = A(R,L m ) + ^(R,L m )p. (20) 

The left hand side of Eq. (|2"0"T) is first order in 1/c 2 . 
Therefore, we will estimate the right hand side of Eq. ([20)) 
in the same order of approximation. By using a Tay- 
lor series expansion to first order of approximation we 
obtain first for the function A (R, L m ) the approximate 
representation 



where 



A{R,L m ) = a + j3R + jL T 



/(0,0) 



and 



6/fi(0,0)' 

fl (0, 0) + /(0, 0)^(0, 0)/ fl (0, 0) 
' 6/3(0,0) 

/HL ra (0,0)/(0,0) 
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6/1(0,0) 



(21) 



(22) 



(23) 



(24) 



respectively. 

For the function (2/3)$ (R,L m ) p we obtain 

2 2 
-$(i?,L m )p« -Sp, 



where 



5 = 



/L m (0,0) 



/«(o,o) ■ 

With these approximations Eq. (|2"0)l becomes 



(25) 



(26) 



(27) 



By substituting into Eq. ([27]) the weak field approxi- 
mation of R, R ps 2A</), gives the generalized Poisson 
equation in modified theories of gravity with an arbitrary 
mattcr-gcomctry coupling as 



a + (3R+ [j+ -5) p. 



. , 7 + 2(5/3 

A0 ps — p + A , 

^ 1 - 2/3 ^ ' 



(28) 



where the condition f3 ^ 1/2 must hold for all p, and 
where we have denoted Aq = a/ (1 — 2/3). 



In order to obtain the correct limit of the Newtonian 
Poisson equation, the function / and its derivatives es- 
timated at the point (0, 0) must satisfy the condition 
(7 + 25/3) / (1 - 2/?) = 1/2. The constant A plays the 
role of an effective cosmological constant, which is natu- 
rally generated in the present model. In most of the as- 
trophysical applications Ao can be neglected. Therefore 
in the first approximation for the potential of the grav- 
itational field of a single particle of mass m we obtain 
cj) (r) = — m/87rr, which is the expression of the standard 
Newtonian potential. As a result, in the equation of mo- 
tion of the test particles one can take a = — V<^ = ajv, 
where a at is the Newtonian acceleration of the particle. 



IV. GEODESIC DEVIATION AND THE 
RAYCHAUDHURY EQUATION WITH AN 
ARBITRARY CURVATURE-MATTER 
COUPLING 

As one can see from Eq. (fTU|) . the proper acceleration 
cPx^/ds 2 is not a covariant object. In particular, its 
vanishing or non-vanishing has no observer-independent 
meaning. In contrast, the relative acceleration between 
worldlines is a covariant quantity, and its vanishing or 
non- vanishing, does not depend on the frame of reference. 

Consider a one-parameter congruence of curves 
x^ (s; A), so that for each A = Ao = constant, x^ (s, Ao) 
satisfies Eq. (fT0|) . We suppose the parametrization to 
be smooth, and hence we can introduce the tangent 
vector fields along the trajectories of the particles as 
= dx^ (s; A) /ds and n M = dx^ (s; A) /dX, respec- 
tively. We also introduce the four-vector 



dx* (s; A) 



dX 



SX = n^SX, 



(29) 



joining points on infinitely close geodesies, corresponding 
to parameter values A and A + <5A, which have the same 
value of s [l^, [2l| . From the definition of and it 
follows that they satisfy the relation dU^/dX = dn^/ds. 
Then it can be easily shown that n v V v U^ = U u V u n t1, 
fig. [2li . Now consider the second derivative 



U V V V (E/ Q V a n") = U V V V {n a V a U^) 
(V„V Q l^) n a U v + (V„n Q ) (V Q C/^) U v . 



D 2 n^ 
~ds^ 



(30) 

By changing the order of covariant differentiation by 
using the definition of the Riemann curvature tensor 

(V r Va-v„v„)tf" = -R% ua ue H3, E q . m 

can be written as 

= K aP n a U p U v +V a (U»V„Unn a . (31) 

By taking into account Eq. (|TU|) . after multiplication 
with the constant factor SX, we obtain the geodesic de- 
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viation equation (Jacobi equation) as [19j, 121142 



D 2 rf 



(32) 



In the case / M = we reobtain the standard Jacobi 
equation, corresponding to the geodesic motion of test 
particles. The interest in the deviation vector derives 
from the fact that if Xq(s) = x^ (s; Xq) is a solution of 
Eq. (|32|) . then to first order Xj"(s) = Xq(s) + •rf' is a solu- 
tion as well, since x^ (s; Ai) ps x^ (s; Ao) +n M (s; A ) <5A pa 
x^ (s;X + 6X). 

Note that the geodesic deviation equation Eq. (f3"2"j) can 
also be derived from the Lagrangian [2l[ 



1 Dif 

L{ll) = 2 9 ^^T^- 



Ir 



(33) 



In this Lagrangian the metric, connection and curva- 
ture are those of a given reference geodesic Xq(s), with 
t/ M (s) = Xq (s) representing the four velocity along the 



same geodesic. These quantities are the background vari- 
ables. The "rf (s) are the independent Lagrangian gener- 
alized coordinates, which are to be varied in the action 
according to the Lagrange equations, 



d 



dL 



dL 



ds d (Dr] v /ds) drf 



(34) 



Then these Lagrange equations give again the equation 
of the geodesic deviation. 

By taking into account the explicit form of the extra- 
force given by Eq. in modified theories of gravity 
with a curvature-matter coupling, the geodesic deviation 
equation can be written as 



+V a V a < V v In 



/i m (R, L m ) 



ds 2 
dL m (p) 



dp 



Explicitly, the geodesic deviation equation becomes 



dL rn (p) 



/i m (R, L m ) 



dL m (p) 



dp 



-rfV^ln 



/i„ (R, L m ) 



dp 



(35) 



As a specific example, consider the case of a linear coupling between curvature and matter [Tl| . where the Lagrangian 
given by f(R, L m ) = fi(R) + Xf 2 (R)L m , and fi(R) (with i = 1, 2) are arbitrary functions of the Ricci scalar R and 
the strength of the interaction between f2{R) and the matter Lagrangian is characterized by a coupling constant . 
For this case, the geodesic deviation equation can be written as 



J»V _ ^ aTT p TT v , A 1 + ) <MR)}d[^F 2 (R)}/dR-l 
ds 2 



R^ru p u v +x 

XF 2 (R) 



[1 + Xf 2 (R)Y 



(L m cr -T^)ifV a RV v R 



1 + A/ a (R) 



{L m gT -T^)ifV a V v R 



XF 2 (R) 
1 + A/ a (R) 



(36) 



Note that a second order tensor Vj,[/ M can be de- 
composed into symmetric and antisymmetric parts, and 
the symmetric part can be further decomposed into a 
trace and trace-free part. Thus, in general, we can write 

V^t/j, = -Oh^ + a fiv + Ufa, + U^Uy, (37) 

where = g^ - UfJJ v , % = U V V V U^ 9 = V ' U U V is 
the expansion of the congruence of particles. The shear 
o^v is given by 

o> = v (^) _ \ 6h ^ - U^U u) , (38) 



where 

V( ll U v) = ^(V„Uf l + V fl U u ), (39) 
and the vorticity uj^ is defined as 

= V { f,U v] -U { fJJ v] , (40) 

respectively, where 

V\* u »] = \ ( v -^ - V^) ■ (41) 

The term V ' fJJ v takes into account the possible presence 
of other forces, which are orthogonal to the four- velocity, 
with four-acceleration given by = U U V ' v l) M . 
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From the definition of the Ricmann curvature we have 

(V„V a - V Q V„) U" = -R$ ua U p . (42) 

By contractin g wi th p = v and after multiplication with 
U a we obtain 

U a V„V a U v - U a V a 6 = -R aP U a U^. (43) 

The first term in this equation can be written as 

U a V v V a U v = V„ (U a V a U v ) - (y v U a ) (V a U u ) . (44) 

Hence we obtain the Rachaudhury equation in the pres- 
ence of an extra force as [23, |23| 



1 



(45) 



where a 2 = o^a^ and oj 2 = uj^oj^, respectively. 

With the use of the field equation Eq. ([6]) , and the ex- 
pression of the extra force, in modified theories of grav- 
ity theories with an arbitrary coupling between curvature 
and matter, the Raychaudhury equation assumes the fol- 
lowing generalised form 
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(a 2 - uj 2 ) + A (R, L, 



V„ < V„ In 



1 



/i„ (R, L m ) 



dL m (p) 



dp 



(L m g" v - T» v ) 



fit (R, L m ) 



U^U u V^ufR(R,L m ) 



1, 



+$(R,L m ) [T^U" - -T 



The vorticity satisfies the equation [22| . \2c 



(46) 



(47) 



while the dynamics of the shear <7 M „ is described by the 
equation [22l. l23l| 

1 2 



u a u A , 



(48) 



where 



C^avX — R)i<jv\+ 9n[\Rv]<j + 9<j[vR\]n H • , (49) 

is the Weyl tensor. With the use of the field equations 
Eq. ([6]) and of the expression of the extra force / M , the 
evolution equations for the shear and vorticity can also be 
obtained explicitly for modified gravity with an arbitrary 
curvature-matter coupling. 



V. TIDAL FORCES WITH AN ARBITRARY 
CURVATURE-MATTER COUPLING 

Tides are the manifestation of a gradient of the gravi- 
tational force field induced by a mass above an extended 
body or a system of particles. In the Solar System tidal 
perturbations act on compact bodies such as planets, 
moons and comets. On larger scales than the solar sys- 
tem, as in a galactic or cosmological context, one can ob- 
serve tidal deformations or disruptions of a stellar cluster 
by a galaxy, or in galaxy encounters 24 \ . In the relativis- 
ts theories of gravitation, as well as in Newtonian grav- 
ity, a local system of coordinates can be chosen, which 
is inertial except for the presence of the tidal forces. In 
strong gravitational fields, relativistic tidal effects can 
lead to interesting phenomena, such as the emission of 
tidal gravitational waves [24| • Relativistic corrections to 
the Newtonian tidal accelerations caused by a massive 
rotating source, such as, for example, the Earth, could 
be determined experimentally, at least in principle, thus 
leading to the possibility of testing relativistic theories of 
gravitation by measuring such effects in a laboratory. 



A. There are no geodesic reference frames in 
f(R,L m ) gravity 

If the Christoffel symbols r™^ associated to a given 
metric are symmetric, it is always possible to chose a co- 
ordinate system in which all the T" become zero at a 
previously assigned point [19J. The corresponding coor- 
dinates are called Riemann normal coordinates, and the 
system of reference can be called a locally Minkowski 
system. Let the given point be the origin of the coordi- 
nate system, and let the initial values of the Christoffel 
symbols be equal to (r^ /3 )o- In the neighborhood of the 
origin we introduce the coordinate transformation 



a/" = x» + I ( rjt „) x a x 



2 \ <*PJo 



(50) 



which provides 



dafdxp){dx f a)~^)o' (51) 



Then, from the transformation law of the Christoffel sym- 
bols, 



r 



lv ( dx^ \ / dx' x \ ( dx 



Xa \dx' v J \dx a ) \dxP 
d 2 x' a \ fdx^\ 



dx a dx@ ) \dx'<T J 



(52) 



it follows that all the transformed Christoffel symbols 



r'£ a are zero. 



One can also show that by a suitable choice of the co- 
ordinate system one can make all the go to zero not 
only in a point, but all along a given world line fl9j |. At 
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the same time with the vanishing of the Christoffel sym- 
bols, the metric g^ v can be reduced, at a given point, to 
its diagonal (Minkowskian) form, g^ v = 77^, and conse- 
quently in the given point the covariant derivatives coin- 
cide with the ordinary partial derivatives, V M = d/dx^ 

In such a coordinate system the equation of motion 
Eq. (flD) takes the form 



d 2 x^ 



dt 2 



n 



(53) 



geod 



where for small particle velocities ds w dt, where t is the 
time coordinate. 

Therefore, in theories of gravity with an arbitrary 
curvature-matter coupling, despite the fact that one can 
locally cancel the Christoffel symbols, and can introduce 
a flat (Minkowskian) metric, due to the presence of the 
extra-force / M , which generally cannot be reduced to 
zero, the motion of the test particles is non-inertial, and 
they will always experience a supplementary acceleration 
induced by the presence of the coupling between matter 
and curvature. 



B. Tidal forces in f(R,L m ) gravity 

In the following we will denote a reference frame in 
which all the Christoffel symbols vanish by a prime. In 
such a system one can always take 77' = 0, which means 
that the particle accelerations are compared at equal 



times. rj n is then the displacement of the particle from 
the origin [25j. Moreover, in the static/stationary case, 
in which the metric, the Ricci scalar and the thermody- 
namic parameters of the matter do not depend on time, 
f'° = 0. With the use of the equation of motion this con- 
dition implies U'° = constant = 1. Therefore, with these 
assumptions, the equation of the geodesic deviation (the 
Jacobi equation) takes the form 



R'owv" + R'ji m v' l U' j U ™ + (54) 



d 2 if 
Ittf 2 

Equation (|54|) can be reformulated as 



F 



dt' 2 jV 1 



(55) 



where F' % is the tidal force, and we have introduced the 
generalized tidal matrix K\ [24j], which is defined as 



QfH 



(56) 



The tidal force has the property dF' % /drj' 3 = Kj, and its 
divergence is given by dF' l /dr]' 1 = K, where the trace K 
of the tidal matrix is 



K = K} = R 00 + R' km U ,k u' m 



dx'i 



(57) 



With the use of the gravitational field equations Eq. 
we can express K as 



K 



1 

+ .fn(R',L'J 
d f d 



1 



2 



f R (R',Ljdt' 2 
d 2 



f R (R\ VJ + 0> (#, L' m ) (t^ - \T'm^j + A (R', L'J Vkm U ,k U' 



jj/kjj/m 



dx' k dx' m 



dx' k I dx 



In 



JL m (R ,L m ) — 



f R {R', L'J + {&, L'J [r' km - l -T m ^j U' k U' m 

I T 1 rJwn rplkrr 



(58) 



Then, since in the Newtonian limit one can omit 
all time derivatives, we obtain for -Rqzo the expression 
R l Ma = d 2 (j) I 'dx 1 dx 1 [25[, where, for simplicity, in the 
following we will omit the primes for the geometrical 
and physical quantities in the Newtonian approximation. 
In Newtonian gravity tu = —d 2 <fi/dx l dx l represents the 
Newtonian tidal tensor (25[. Therefore, in modified grav- 
ity with a curvature-matter coupling wc obtain the tidal 
acceleration of the test particles as 



dt 2 



= F l 



if- 



dx l dx l 



dx 1 



(59) 



where V 3 and V m are the Newtonian three-dimensional 



velocities. In the Newtonian approximation, in modified 
theories of gravity with a curvature-matter coupling the 
tidal force tensor is defined as 



and its trace gives the generalized Poisson equation, 



w = A* + R jm VW- + J-. (61) 
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C. The Roche limit in modified gravity with an 
arbitrary curvature-matter coupling 

In Newtonian gravity, the spherical potential of a given 
particle with mass M is cf>(r) = —M/8wr. By choosing a 
frame of reference so that the x-axis passes through the 
particle's position, corresponding to the radial spherical 
coordinate, that is, (x = r,y = 0,z = 0), the Newto- 
nian tidal tensor is diagonal, and has the only non-zero 
components 



= diag 



2M 
8nr 3 



M 



M 



■ r 3 l 



(62) 



The Newtonian tidal force F t can be written as 



F tx = 2MAx/87rr 3 , F t 



-MA; 



ty 



-GMAy/8irr 3 and F u 



respectively [25[. These results can be 



used to derive the generalization of the Roche limit in 
modified gravity with an arbitrary coupling between mat- 
ter and curvature. 

The Roche limit is the closest distance Troche that a 
celestial object with mass m, radius R m and density p m , 
held together only by its own gravity, can come to a mas- 
sive body of mass M, radius Rm and density pm, respec- 
tively, without being pulled apart by the massive object's 
tidal (gravitational) force [lij . For simplicity we will con- 
sider M ^ to, so that the center of mass of the system 
nearly coincides with the geometrical center of the mass 
M. 

The elementary Newtonian theory of this process is as 
follows. Consider a small mass Am located at the surface 
of the small object of mass m. There are two forces acting 
on Ato, the gravitational attraction of the mass m, given 
by 



Fc 



toAto 

87Ti?2 ' 



(63) 



and the tidal force exerted by the massive object, which 
is given by 



Ft 



MAmR 
8irr 3 



(64) 



where r is the distance between the centers of the two ce- 
lestial bodies. The Roche limit is reached at the distance 
r = rR OC h ei when the gravitational force and the tidal 
force exactly balance each other, Fq = Ft, thus giving 

241 



1/3 / \l/3 

- ) = 2 1 ' 3 R M i^-J . (65) 



In modified gravity with a curvature-matter coupling 
the equilibrium between gravitational and tidal forces oc- 
curs at a distance rR OC h e given by the equation 



ft I 



, R r jrm ViV m + ?£) Rm = + r, 

*™ Roche dr J 87ri? ". 



where f r is the radial component of the extra-force, 
which modifies the Newtonian gravitational force, and 
the curvature tensor R^ rm (no summation upon r) must 
be evaluated in the coordinate system in which the 
Newtonian tidal tensor is diagonal. Hence we obtain 
the generalized Roche limit in the presence of arbitrary 
geometry-matter coupling as 



M 

T Roche ~ Rm I 

TO 



1/3 



8nR 3 

3to 



df r 

3™ 1 g r 



" | 87T 

3m 



(67) 



(66) 



where we have assumed that the gravitational effects due 
to the coupling between matter and curvature are small 
as compared to the Newtonian ones. 



VI. CONCLUSIONS 

In the present paper we have developed some of the ba- 
sic theoretical tools necessary to investigate the proper- 
ties of gravitational models with direct curvature-matter 
interaction. In particular, we have obtained the Ray- 
chaudhury equation, which can be used to analyze the 
existence of black holes in these types of theories, and to 
rigorously formulate the singularity theorems. The pres- 
ence of the extra-force in the gravitational model leads to 
the appearance of some new terms in the Raychaudhury 
and the geodesic deviation equation, terms which are 
a direct consequence of the curvature-matter coupling. 
Therefore the presence of the extra force can have a sig- 
nificant effect on the formation of massive astrophysical 
objects by gravitational collapse. However, we have to 
mention in the present gravity model the full non-linear 
theory contains higher than two spacetime derivatives in 
the field equations. Hence modified gravity theories with 
geometry-matter coupling can be considered as effective 
theories only apart from the special Einstein case. 

The geodesic equation can be used to study the effects 
of the generalized tidal forces, which could lead to the 
possibility of observationally testing the model through 
the observational effects of the tides due to an extended 
mass distribution. Typical situations in which the effects 
of the tides are of major importance are galactic encoun- 
ters, globular clusters under the influence of the galactic 
mass distribution, and Oort cloud perturbations by the 
galactic field (26|. It was suggested that the compres- 
sive tidal field in the center of the fiat core of early type 
galaxies and of the ultra-luminous galaxies compresses 
the molecular clouds, producing the dense gas observed 
in the center of these galaxies. The curvature-matter in- 
teraction modifies the nature of the tidal forces. There- 
fore the observational study of the tidal forces may give 
some insights in the fundamental aspects of the gravita- 
tional interaction. 
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On the other hand the Newtonian limit of the geodesic 
deviation equation Eq. (|59[) shows that the tidal motion 
of test particles is directly influenced not only by the gra- 
dient of the extra force /, which is basically determined 
by the gradient of the Ricci scalar, but also by an explicit 
coupling between the velocity and the Riemann curva- 
ture tensor. The possibility that the Pioneer anomaly 
[l?} is due to the presence of the extra force _due to the 
geometry-matter coupling was analyzed in 



111. How- 



ever, Eq. ([5^| shows that this effect may also be velocity 
dependent, leading to the possibility of the detection of 
these type of extra gravitational effects by the method of 
the Doppler tracking of a spacecraft. On the other hand, 
recent researches on the Pioneer anomaly |28| have given 



convincing arguments against its modified gravity origin. 

To conclude, the present study opens further possibili- 
ties for the experimental and observational investigation 
of the possible coupling between curvature and matter, 
and for a large class of modified gravity theories. 
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